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The NA48/2 experiment at the CERN SPS has collected an unprecedented sample of — ► 37T 
decays. The high statistics and the good resolution of the detectors allow a unique investigation of 
the detailed phase space distributions of these decays. The effects of final state pion rescattering 
observed in the Dalitz plot distribution of the — » n ± n Q n Q decays turned out to be a powerful 
tool for extraction of the S-wave nn scattering lengths. The large statistics also allowed a precise 
measurement of the Dalitz plot slope parameters for the K — > 3n ± decays. 
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Introduction 

The primary goal of the NA48/2 experiment at the CERN SPS is the search for direct CP 
violation in — > 371 decays [jl]]. Data have been collected in 2003-04, providing samples of 
~ 4 x 10 9 fully reconstructed — > 37T 111 and ~ 10 8 ^T* — > tt*?! ?! decays. Surprisingly, a study 
of a partial sample of — > ^n ^ decays corresponding to about 25% of the total sample 
revealed an anomaly in the 7l°n° invariant mass (Moo) distribution in the region around Moo = 2m + , 
where m + is the charged pion mass [0|. This anomaly, dubbed "cusp effect", never observed in 
previous experiments, was theoretically interpreted as an effect due mainly to the final state charge 
exchange scattering process 7r + 7T~ — ► Tt % in — > 3n ± decay, and was shown to provide a 
precise determination of — a%, the difference between the 5- wave nn scattering lengths in the 
isospin 1 = and 1 = 2 states [J3]] . A number of theoretical approaches to describe this process are 
being developed; the original NA48/2 measurement of ao — «2 was performed in the framework of 
the approach [Q]. The current paper presents a new step of the analysis, namely a preliminary result 
of a measurement based on the full NA48/2 data sample within the same theoretical framework. 

In addition, a measurement of the Dalitz plot slopes of the — > 3%^ decay based on a partial 
NA48/2 data sample is presented. 

1. NA48/2 experimental setup 

Two simultaneous K + and K beams are produced by 400 GeV protons impinging on a 40 cm 
long Be target. Particles with a central momentum of 60 GeV/c and a momentum band of ±3.8% 
produced at zero angle are selected by a system of dipole magnets forming an "achromat" with null 
total deflection, focusing quadrupoles, muon sweepers and collimators. With 7 x 10 11 protons per 
burst of 4.5 s duration incident on the target the positive (negative) beam flux at the entrance of the 
decay volume is 3.8 x 10 7 (2.5 x 10 7 ) particles per pulse, of which 5.7% (4.9%) are K + (K~). The 
decay volume is a 1 14 m long vacuum tank. 

Charged particles from decays are measured by a magnetic spectrometer consisting of four 
drift chambers and a large-aperture dipole magnet located between the second and third chamber. 
Each chamber has eight planes of sense wires: two horizontal, two vertical and two along each 
of two orthogonal 45° directions. The spectrometer is located in a tank filled with helium at at- 
mospheric pressure and separated from the decay volume by a thin (0.31%Xo) Kevlar window. A 
16 cm diameter vacuum tube centered on the beam axis runs through the spectrometer and subse- 
quent detectors. Charged particles are magnetically deflected in the horizontal plane by an angle 
corresponding to a transverse momentum kick of 120 MeV/c. The momentum resolution of the 
spectrometer is a{p)/p = 1.02% © 0.044% p (p in GeV/c). The spectrometer is followed by a 
scintillator hodoscope consisting of two planes segmented into horizontal and vertical strips. 

A liquid krypton calorimeter is used to reconstruct — ► n ± 7l Q 7i (> decays. It is an almost 
homogeneous ionization chamber with an active volume of 7 m 3 of liquid krypton, segmented 
transversally into 13248 projective cells of 2x2 cm 2 by a system of Cu-Be ribbon electrodes, and 
with no longitudinal segmentation. The calorimeter is 27Xo thick and has an energy resolution 
o{E)/E = 0.032//E©0.09/£© 0.0042 (E in GeV). Spatial resolution for a single electromag- 
netic shower is a x = a y = 0.42%/\/Z?©0.06 cm for each transverse coordinate x,y. 
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Figure 1: Reconstructed spectra of n°7i° invariant mass showing evidence for the cusp effect: the full 
kinematic ranee for (a) 2003 data (16.0 x 10 6 events), (b) 2004 data (43.6 x 10 6 events); zoomed threshold 
region for (c) 2003 data, (d) 2004 data. The 2003 plots correspond to the original discovery of the effect [g], 
while the 2004 plots correspond to the progress with respect to the original analysis. 

A detailed description of the components of the NA48 detector can be found elsewhere [||]. 

2. Cusp effect and measurement of pion scattering lengths 

The reconstructed spectra of 71° 71° invariant mass Moo for 2003 and 2004 data samples (totally 
59.6 x 10 9 events) are presented in Fig. [[]. The change of slope at 7t + n~ threshold is clearly visible. 
For description of this effect the K^ 1 — > 71*71° 71° amplitude is presented as a sum of two terms: 



''1) 



(2.1) 



where ^#o is the "unperturbed" amplitude expressed as a polynomial expansion in terms of the 
kinematic variables u = (53 — so)/m 2 + and v = (s\ — S2)jnv\, where = (Pk — Pi) 2 , sq = (si +S2 + 
*3)/3, Pk and Pi are 4-momenta of kaon and pions, and i= 1,2 correspond to the two "even" (i.e. 
identical) pions: 

(2.2) 



?o(u,v) = ^o(0,0) • (1 +g u/2 + h'u 2 /2+kV/2) 



U..2 
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Figure 2: Deviation of the data spectrum from the fit result with statistical errors (combined 2003+2004 
data set): A = Data/ Fit — 1. Good quality of the fit and an excess of events in the region of the threshold are 
demonstrated. 

and is a contribution from the — > 3^ decay amplitude through 7i + 7i~ — > 7r°7r° charge 
exchange, which in particular simplest case of the original Cabibbo theory [^] is given by 

Jt\ = -2a x m + J? + J 1 - (M 00 /2m + ) 2 . (2.3) 

Here, in the limit of exact isospin symmetry, a x = («o — fl2)/3. The amplitude ^\ changes from 
real to imaginary at the threshold Moo = 2m + ; as a consequence it interferes destructively with ^#o 
below the threshold (leading to 13% integral depletion in this region), and adds quadratically above 
the threshold. 

The model used for the present measurement is based on the formulation [Q], which takes 
into account all rescattering precesses at the one-loop and two-loop level. In this approach the 
matrix element of the — > n ± n Q 7i° decay includes a number of additional terms depending on 
five 5-wave %% scattering lengths (corresponding to the processes n + 7i~ — ► 7r°7r°, n + 7i + — ► n + 7i + , 
— > 7r + 7T _ , 7T + 7T —> 7i + 7t° and 71° 71° — ► 71° 71°) expressed as linear combinations of ao and <32- 
In addition to [Q], isospin breaking effects are taken into account introducing a single parameter 
£ = (rn 2 + - ml) j m \ = 0.065 [§]. 

The fit to extract the scattering lengths and Dalitz plot slopes go, h! was performed in the Moo 
projection of the data using a full GEANT-based Monte Carlo simulation of the detector response. 
The used rescattering model does not include radiative corrections, which are particularly important 
at the threshold Moo = 2m + , and contribute to formation of 7r + 7T~ atoms (pionium). Thus a group 
of seven bins near the threshold has been excluded from the fit. The quality of the fit (# 2 /NDF = 
164/139 for 2003 analysis, and X 2 /NDF = 119/139 for 2004 analysis) illustrated in Fig. | shows 
an excess of events in this excluded region. This excess, being interpreted as due to pionium 
formation, yields the rate of pionium formation R = T(^ ± — ► 7T + A2^)/r(A r± — >37T ± ) = (1.82± 
0.21) x 10~ 5 , somewhat higher than a theoretical prediction 0. 

Measurement of the quadratic Dalitz plot slope k! was performed using the v projection of the 
data and fixing the values of ao, ci2, go and h' measured by the above method. Then the fit in Moo 
projection was re-iterated to account for the measured non-zero value of k! . 
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Systematic uncertainties due to fitting technique, trigger efficiency, description of geometric 
acceptance and resolution, calorimeter non-linearity, and simulation of showers in the calorimeter 
have been evaluated. External uncertainties due to limited experimental knowledge of / 
at the % + K~ threshold have been also considered. Stability checks with respect to decay vertex 
position, particle separations in the calorimeter front plane, and kaon sign have been performed. 

Results and conclusions 

The original NA48/2 measurement of the of %k scattering lengths [||] by exploring the cusp 
effect in the K — > t^tt ?! decay spectrum has been improved: the full NA48/2 data sample has 
been used, and a more elaborate study of systematic effects performed. The model [Q] with isospin 
breaking corrections has been used. The measured scattering lengths are: 

(a -a 2 )m + = 0.261 ± 0.006 stat . ± 0.003 syst . ± 0.001 ex t. 
a 2 m+ = -0.037 ± 0.013 stat . ± 0.009 syst . ± 0.002 ext . 

The external uncertainties are due to the limited knowledge of r(/T ± — > n ± 7l (> 7i Q ) /^(K^ — > 371*). 
Moreover, an uncertainty m + d(ao — a 2 ) = 0.013 has to be attributed to the result due to precision 
of the theoretical model. The Dalitz plot slopes corresponding to the used model are found to be 

go = 0.649 ± 0.003 stat . ± 0.004 syst . 
h' = -0.047 ± 0.007 stat . ± 0.005 syst . 
kf = -0.0097 ± 0.0003 stat . ± 0.0008 syst . 

In addition, a measurement of the Dalitz plot slopes of the PDG parameterization [Jsj] of the — > 
371^ decay with a sample of 4.71 x 10 8 fully reconstructed events yielded the following results: 

g = -0.21134 ± 0.00013 stat . ± 0.00010 syst . 
h = 0.01848 ± 0.00022 stat . ± 0.00033 syst . 
k = -0.00463 ± 0.00007 stat . ± 0.00012 syst . 

This measurement is described in detail in [^]. The results are compatible with the world average, 
and demonstrate the validity of the conventional parameterization at the new level of precision. 
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